NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



REPORT No. 665 

CALCULATION OF THE AERODYNAMIC 
CHARACTERISTICS OF TAPERED WINGS WITH 
PARTIAL-SPAN FLAPS 

By HENRY A, PEARSON and RAYMOND F. ANDERSON 



I 




193 9 



REPRODUCED BY I 

NATIONAL TECHNICAL 

INFORMATION SERVICE | 0 

U.S. DEPARTMENT OF COMMERCE > 
SPRINGFIELD, VA. 22161 j 



J " 



AERONAUTIC SYMBOLS 



1. FUNDAMENTAL AND DERIVED UNITS 





Symbol 


Metric 


English 


. Unit 


Abbrevia- 
tion 


Unit 


Abbrevia- 
tion 


Length 


I 
t 

F 


meter 


m 

s 

k S 


foot (or mile) ___ 


ft. (or mi.) 
sec. (or hr.) 
lb. 


Time____ 

Force 


second _ _ 

weight of 1 kilogram _ _ 


second (or hour). 

weight of 1 pound. 


Power, . 


P 

V 


horsepower (metric) 




horsepower. __ 


hp. 

m.p.h. 

xf.p.s. 


Speed. 


f kilometers per hour 

[meters per second 


k.p.h. 
m.p.s. 


miles per hour. . 

feet per second___i_--_ 



2. GENERAL SYMBOLS 



W, 

g, 

I, 



Weight =mg 

Standard acceleration of gravity =9. 80665 
m/s 2 or 32.1740 ft./sec. 2 

■KM ■ W 

Mass— — 
9 

Moment of inertia— mk 2 . (Indicate axis of 

radius of gyration k by proper subscript.) 
Coefficient of viscosity 



v } Kinematic viscosity 

p, Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-nr 4 -s 2 at 

15° C. and 760 mm; or 0.002378 'lb.-ft .-* sec. 2 
Specific weight of "standard" air, 1.2255 kg/m 3 or 

0.07651 lb./cu. ft. 



3. AERODYNAMIC SYMBOLS 



s, 
o, 

h, 

c, 
P 

s' 
v, 

a, 
d, 

Do, 
D i} 

D P: 

C, 
R, 



Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio 

True air speed 

Dynamic pressure =^pV 2 

L 

Lift, absolute coefficient (7^= —g 
Drag, absolute coefficient C D —-^L 



D 0 



Profile drag, absolute coefficient @r> 0 —-^g 



Dt 



Induced drag, absolute coefficient G Di =-^ 



D r 



Parasite drag, absolute coefficient G Dp =^-^ 



Cross-wind force, absolute coefficient (7 C = 
Resultant force 



VI 
p — > 



i m Angle of setting of wings (relative to thrust 
line) 

i t , Angle of stabilizer setting (relative to thrust 
line) 

Q, Resultant moment 

ft, Resultant angular velocity 

Reynolds Number, where I is a linear dimension 
(e.g., for a model airfoil 3 in. chord, 100 
m.p.h. normal pressure at 15° C, the cor- 
responding number is 234,000; or for a model 
of 10 cm chord, 40 m.p.s., the corresponding 
number is 274,000) 
Center-of-pressure coefficient (ratio of distance 

of c.p. from leading edge to chord length) 
Angle of attack 
Angle of downwash 
Angle of attack, infinite aspect ratio 
Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 
7, Flight-path angle 



Gn 



<*0, 
<X a , 



/ 
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CALCULATION OF THE AERODYNAMIC CHARACTERISTICS OF TAPERED 

WINGS WITH PARTIAL-SPAN FLAPS 

By Henry A. Pearson and Raymond F. Anderson 



SUMMARY 

Factors derived from wing theory are presented. By 
means of these factors, the angle of zero lift, the lift-curve 
slope, the pitching moment, the aerodynamic-center posi- 
tion, and the induced drag of tapered wings with partial- 
span flaps may be calculated. The factors are given for 
wings of aspect ratios 6 and 10, of taper ratios from 0.25 
to 1.00, and with flaps of various lengths. 

An example is presented of the method of application of 
the factors. Fair agreement with experimental results is 
shown for two wings of different taper ratio having plain 
flaps of various spans. 

INTRODUCTION 

Because of the widespread use of tapered wings 
equipped with partial-span flaps, it is desirable to have 
means for computing their aerodynamic characteristics. 
Previous reports (references 1, 2, and 3) have presented 
theoretical factors for use in computing the aerodynamic 
characteristics of wings with linear and with arbitrary 
twist and for use in finding the load distribution of 
wings with partial-span flaps. 

This report presents factors, based on airfoil theory, 
for use in calculating the induced drag, the angle of 
zero lift, the pitching moment, and the aerodynamic 
center of tapered wings with partial-span flaps of con- 
stant flap-chord ratio. The factors, when used with 
adequate section data, should apply to various types of 
flap and various amounts of flap deflection. 

THEORETICAL RESULTS 

The particular wing chord distributions for which the 
theoretical computations were specifically made are 
given in figure 1 where the wing quarter-chord line is 
shown as straight. Two aspect ratios (A=6 and 10) 
and four taper ratios (X = 1.00, 0.75, 0.50, and 0.25) 
were used. A list of the symbols used herein is given 
in appendix A. Inasmuch as the various characteris- 
tics for elliptical wings with partial-span flaps could be 
obtained relatively easily, they were sometimes com- 
puted in order to aid in determining the shape of the 
various computed curves for the tapered wings. 

The span load distributions from which the aero- 
dynamic characteristics were obtained are given in ref- 




'M 

25 




Figure 1.— Wing chord distributions. 



erence 2 where a slope of the section lift curve equal to 
5.67 per radian was used. The computations apply 
only to those cases in which no aerodynamic twist is 
present before the flaps are deflected. 

Although the ordinary lifting-line theory is applica- 
ble only to wings without sweepback, experimental evi- 
dence indicates that small amounts of sweepback 
have no appreciable effect on the span loading. The 
computations may thus be applied to wings with moder- 
ate sweepback as long as the chord distributions are 
similar to those indicated in figure 1 . 

The computed aerodynamic characteristics are given 
in terms of factors such as J, H, and G. The method of 
calculating the factors is omitted because of its length, 
but the formulas for the factors are presented in ap- 
pendix B. The physical significance of the factors and 
of the aerodynamic characteristics they represent, how- 
ever, is explained in the following sections. 

Angle of zero lift. — The change in the angle of zero 
lift of a finite wing accompanying a flap deflection de- 
pends upon several variables, such as flap span, flap 
deflection, flap chord, and flap type. The effect of the 
last three variables can be conveniently represented 
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by the section characteristic Ac h the increment of 
section lift coefficient obtained by deflecting the flap. 

By this grouping of variables, the change in the angle 
of zero lift (in degrees) for a wing with partial-span 
flaps can be expressed by the equation 

Aa S(£ta0) =— JAci (1) 

In order to obtain the angle of zero lift for the wing, 
this increment must be added to the initial angle of 
zero lift, i. e., the angle before the flap is deflected. 
If this initial angle is measured from the chord of the 
root section, as is usually the case, the angle of zero 
lift for the wing is given by 

■JA Cl (2) 




.2 .4 .6 

Flap span/wing span, 



b f /b 



FjCiUUE 2.— Factor of angle of zero lift, /. 



The computed variation of the factor J with flap span 
is shown in figure 2 for various aspect ratios and taper 
ratios. 

Although the values of J given in figure 2 apply 
specifically to wings in which the flap-chord ratio, or 
Ac i, is constant along the portion with flaps and in 
which the flaps begin either at the center or at the tips, 
the results may be used to predict the angle of zero lift 
for any starting point of the flaps and for any Ac t 
distribution as long as they are symmetrical about the 
wing center. For example, if flaps of uniform flap- 
chord ratio extend from 0.36 to 0.76, the proper 
value of the factor J is the difference between the values 
for 0.36 and 0.76 as shown by AJ in figure 2. The 
extension to the case of a nonuniform symmetrical dis- 
tribution of Aci consists simply in considering the 
resulting Ac t distribution to be caused by a series of 
elemental flaps of various lengths and performing either 
a numerical or a graphical integration for the value of J. 
(See procedure given in reference 3.) In cases where 
the variation of Ad along the span is slight, however, 
the use of an average value of Aci is justified. 

lift-curve slope.-— The wing lift-curve slope, a, per 
degree may be found from the equation 



1 + 



57.3a 0 



(3) 



where 

A is wing aspect ratio, b 2 jS 
1.00 



.99 



f 



.98 



.97 



.96 









































A 



























































































































.2 



Figure 3. 



.4 .6 .8 

Taper ratio. A 

-Factor of wing lift-curve slope, /. 
a. 

a=f- 



1.0 



57.3a. 



1+ 



xA 



/, a theoretical factor given in figure 3. This factor 
has been plotted from results given in reference 1. 

~a Q , the weighted average of the section lift-curve 
slopes. 

An average slope, weighted according to chord length, 
must be used because the slope of the sections with flaps 
may be considerably different from the slope of the 
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(4) 



sections without flaps. If the section lift-curve slopes 
are constant across the spans of the flapped and the 
unflapped parts of the wing, then"a 0 may be found in 
terms of the fraction of the area of the wing equipped 
with flaps: 

where 

a 0 is the lift-curve slope of section without flaps, per 
degree. 

a 0f , lift-curve slope of section with flaps, per degree; 
S, area of wing. 

S f , area of part of wing equipped with flaps. 
If a 0 and a 0f are not constant across the two parts of 
the wing, then a 0 may be found by integration. 



namic centers is unaltered by deflecting the flaps, the 
x position of the wing aerodynamic center would be 
the same with the flaps either deflected or neutral. 

If the aerodynamic center of the root section is taken 
as a reference point and the aerodynamic centers of all 
the wing sections are assumed to lie on a straight line 
making an angle A with the lateral axis (see fig. 4), 
then the x location of the load center is given by 

x a . c =Hb tan A (5) 

From equation (5), the aerodynamic-center position 
may also be related to the mean chord S/b and to the 
aspect ratio A by the equation 



S/b 



-HA tan A 



(6) 



Sweepback 
angle A \ >^"" 
\ 

\ 



^^Line of aerodynamic centers 



>" 6f , flap angie 

„ Basic span load distribution 

V , Resultant of upward loads 




Resultant of 
downward loads 



Lateral axis' 

Figure 4.— Typical wing-flap combination showing basic and additional distributions 



"Additional span load 
distribution 



Aerodynamic-center location. — The aerodynamic 
center of a wing is defined as the point about which 
the pitching moment is constant up to high lift coeffi- 
cients; and, since the resultant lift must act through 
this point, the x position of the wing aerodynamic cen- 
ter is, in effect, nothing more than the fore-and-aft 
location of the centroid of the load distribution. 

The load distribution of a wing with flaps is consid- 
ered in this report to consist of the two components 
that are shown in figure 4. The basic load distribution 
is the span loading for zero lift with the flaps deflected; 
its ordinates are proportional to the value of Ac t . The 
additional load distribution is that for the wing with 
flaps neutral; the total lift, however, is the same as 
that for the wing with flaps deflected. Since the basic 
load distribution contributes no lift, it does not enter 
into the determination of the aerodynamic center. 
Thus, if the chordwise position of the section aerody- 



For a wing with sweepforward, A is negative and the 
aerodynamic center of the wing is ahead of the aero- 
dynamic center of the root section. 

Values of H are shown in figure 5. 



.24 



H 



.22 



.20 





















A 




















.10^ 
























£it/p 


tical 




















<£- 


















V 



















.2 .4 .6 .8 

Taper ratio, A 

Figure 5.— Factor of wing aerodynamic center, H. 
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Pitching moment. — As shown in figure 4, the upward 
and the downward parts of the basic load form a couple 
having a magnitude that increases directly with the 
semispan length, the angle of sweepback, and the flap 
deflection, i. e., Ac t . An equation for the pitching 
moment due to the basic load distribution can thus 
be written: 

b 

M ib =lc^Ac l tan A %S 

where k accounts for variations with wing taper, aspect 
ratio, and flap span. Because the basic load distribu- 
tion is zero with no flap and is also zero with a full- 
span flap, the factor k would have a maximum value 
at an intermediate flap span. 

Transforming the preceding equation into the co- 
efficient form gives 

C mib -=GAc l A tn,n A (7) 

where values of G are given in figure 6. For a wing 
with sweepback (A positive), the sign of the pitching- 
moment coefficient due to the basic lift is positive if 
the flap deflection introduces an effective washout 
toward the tip (e. g. , flaps at the center deflected down- 
ward or flaps at the tip deflected upward). For a 
wing with sweepforward, the sign of C mtj) is negative 
for the same flap deflections. 



Flap span/wing span, bf/b (flaps at tip J 
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Flap span/wing span, bf/b (flaps at center) 

Figure 6.— Factor of basic lift pitching moment, G. 
Cm, =GAciA tan A 
l b 

In case the aerodynamic centers do not lie on a 
straight line so that the angle of sweepback is not con- 
stant along the span, C mz& may be graphically obtained 

from the equation 



C-n 



2b 

'S 2 



cm 

Jo 



xc lb cdy 



(8) 



where, at any point along the span, 

x is the moment arm measured from the aerodynamic 
center of the root section and parallel to the root 
chord (positive, rearward; negative, forward). 



Ci b , the section lift coefficient for the basic loading, 
c, the chord. 

In order to obtain the total pitching-moment 
coefficient, that due to the sections must also be added. 
This pitching-moment coefficient is given by 



n 2b P /2 2^ 



(9) 




.4 .6 
Taper ratio, A 

Figure 7.— Factors of section pitching moment, E and E'. 
C n =Ec m .+E'Ac m 



For wings with flaps, however, the value of the section 
pitching-moment coefficient c ma c may be assumed to 
consist of two parts: One denoted by c Wq , the section 
coefficient with flaps neutral; and the other denoted by 
Ac mt the increase in the section coefficient above c m(l 
due to the flaps. If c tllQ is constant across the span and 
Ac m is constant across the flap span (i. e., the flap-chord 
ratio is constant), then the pitching-moment coefficient 
due to the sections can be given by 

C ms =Ec mQ +E'Ac m (10) 

Values of E and E' for these conditions are given in 
figure 7 for the tapered wings. These values have been 
determined from the relations 



rrr 2b r b » 2 , 

77V 2b fW2 2 , 

E= s>L cdy 



If neither c mo nor Ac m were constant across the span, 
then it would be necessary to use equation (9) and to 
evaluate C ms by an integration, as will be illustrated 
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later. The total wing pitching-moment coefficient is 
given by 

^ m a.e.~^ m s Jr ^ m i b (11) 
The coefficient Cm Mm is defined by the equation 



S 2 



(12) 



where M is the total pitching moment. 

Induced drag. — For any wing with a twist that is 
symmetrical about the wing center line, the induced- 
drag coefficient may be given by the equation 



(13) 



The factors u, v, and w for wings with partial-span 
flaps, i. e., for the case of an abrupt twist, are given in 
figure 8. 

The first term on the right-hand side of equation (13) 
is the usual induced-drag coefficient of an untwisted 
wing and the other two terms result from the aerody- 
namic wing twist introduced by deflecting the flaps. 
It can be seen from figure 8 that, for certain taper 
ratios, the v and the w factors are of opposite sign and 
their contributions counteract each other. In fact, under 
certain conditions, the sum of the last two terms may 
be slightly negative; and, as a result, the elliptical 
wing induced-drag coefficient may be approached. 
This tendency exists when the flaps are so placed and 
deflected that an elliptical loading is approximated. 

EXPERIMENTAL RESULTS 

APPARATUS AND TESTS 

In order to provide a check on the reliability of the 
theoretical factors that have been presented, two 
tapered wings with partial-span flaps were tested. In 
addition, tests were made of three rectangular wings 
with full-span flaps to provide section data for use in 
calculating the characteristics of the tapered wings. 
The wings were made of aluminum alloy and had an 
area of 150 square inches. 

A list of the tapered wings and the different flap 
lengths used is given in table I, together with the 
taper ratio, the aspect ratio, and the airfoil sections 
of the root and the construction tip (the extreme tip). 
The tips were rounded as shown in figure 1. The 
N. A. C. A. 23012 tapered wing had a moderate sweep- 
back (line through quarter-chord points) but the N. A. 
C. A. 5-10-16 tapered wing had no sweepback. In the 
construction of the wings, straight-line elements were 
used between corresponding points of the root and the 
construction tip sections. For the N. A. C. A. 23012 
wing, the chords of all sections along the span were in 
one plane; whereas, for the N. A. C. A. 5-10-16 wing, 
the highest points of the upper surface of each section 
were in one plane. The ordinates of the N. A. C. A. 
5-10-16 wing are given in reference 4. 



1.00 



.98 



.96 



.94 



A = 6 
A=IO 



.2 



.6 



1.0 



.4 .6 
Toper ratio, A 

F/ao span I 'wing span, b f /b (flaps at fip) 
1.0 .8 .6 .4 .2 O 



.004 



.002 



§-.002 



-004 



-.006 



-.008 
0 



.008 































— ^ — - 




A = 

j 


25 
















\ 


V \ 
\ \ 




















\^ 

\ 


\ 














Elh 


ptico 


A • 


\ v 














\\ 








"'/£ 

// 

1 


An 










\ 




.50 




-ft— > 
/ 1 /, 

1 /'> 


771 
/i 

i — 








— 






v. 




rfrt 

n 








\ . 


V\ 
V s 
\ 


^.75 








1 1 








^ 




\ 






fl 
// 


~l — 
/ 

/ 










\\ 
\\ 




v. 




1 


/ 






J 

J 


1 = 5 

l = /0 


\ J 

\ 


\ 






T 

/ 

/ 














\ 

x 


WO 


/ 


~? 

i 







-.004 



-.002 



.go?! 



.004 



.006 



.008 



Flap span/wing span, bf/b (flaps at center) 




.2 .4 .6 .8 

FJap span J wing span, bf/b 

Figure 8.— Factors of induced drag, it, v, and w. 
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The rectangular wings had N. A. C. A. 23009, 23012, 
and 23015 sections and were included to provide airfoil 
section characteristics to aid in calculating the charac- 
teristics of the tapered wings. 

Plain 0.2c flaps deflected downward 20° were built 
into all the wings and were made to simulate flaps 
pivoted about the midpoint of the thickness at 0.8c. 
Fillets of small radii were used to join the flap to the 
wing and to seal the gap, as indicated at the top of 
figure 9. 

All the wings were tested in the variable-density 
wind tunnel at a pressure of 20 atmospheres. The lift, 
the drag, and the pitching moment were measured at 



c m(a c j are given about the aerodynamic-center posi- 
tion with the flap neutral. 

The results of the tests of the tapered wings are 
presented in the usual manner in figures 12 to 17. In 
addition to the usual characteristics, the lift-curve 
peaks are given for two values of the effective Reyn- 
olds Number to indicate the scale effect on C Lma ' x . 
The Reynolds Number is based on the mean chord S/b. 
On the right side of the figures, effective profile-drag 
coefficients are given. This coefficient is the total drag 
coefficient with the induced-drag coefficient for ellip- 
tical span loading deducted, that is, Cj> e —C D —C L 2 /wA. 
The values of C D . have been corrected to effective 
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Figure 9.— The N. A. C. A. 23009 airfoil with 0.2c plain flap down 20° 



the usual high Reynolds Number and, in addition, the 
maximum lift was measured at a lower Reynolds Num- 
ber to indicate the scale effect on C Lm . The method 
of making and correcting the tests and a description 
of the tunnel are given in reference 5. 

The results of the tests are presented in the usual 
form as figures 9 to 17. The results of the tests of 
rectangular wings, plotted on the left side of figures 
9 to 11, have been corrected to aspect ratio 6; whereas 
the results given on the right side have been corrected 
to airfoil section characteristics by the method explained 
in reference 6. The pitching-moment coefficients 



Reynolds Number B e by subtracting an increment 
(0.0011) to allow for the reduction in skin-friction drag 
when extrapolating from test to effective Reynolds 
Number (reference 6). The pitching-moment coeffi- 
cients given are based on the mean chord S/b so that 
C m =M/aS(S/b)=Mb/qS 2 . The coefficients for each 
wing-and-flap combination are given about" an axis 
through the aerodynamic center determined by the 
method given in the appendix of reference 4. The 
location of the aerodynamic center given in the upper 
part of the right side of the figures is measured from 
the quarter-chord point of the root chord and is in 
terms of S/b. 
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Where tested: LMA.L.Jest.-V.D.T. 1509 
R e : 8,130.000 Date: 7-/9-37 
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Figure 10.— The N. A. C. A. 23012 airfoil with 0.2c plain flap down 20 
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Ficjuke 11.— The N. A. C. A. 23015 airfoil with 0.2c plain flap down 20 
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Figure 12.— The tapered N\ A. C. A. 230X2 airfoil. 




Figure 13. -The tapered N. A. C. A. 23012 airfoil with 0.36 plain flap down 20°. 
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Figure 18.— Variation of section data with thickness. The N. A. C. A. 230 series airfoils; effective Reynolds Number, 8,200,000. 



CALCULATED CHARACTERISTICS -OF THE WINGS 

The factors previously presented were applied to the 
calculation of the characteristics of the wings used in 
the tests and the results are summarized in table I. 
The calculations will be illustrated for the tapered 
N. A. C. A. 23012 wing with the 0.56 flap. 

Angle of zero lift and lift-curve slope. — The angle of 
zero lift by equation (2) is: 

«'a=o> = -1.2- (6.07 X 0.90) - - 6.7° 

The value of Acj (0.90) was measured from figure 10 at 
approximately the average lift coefficient of the basic 
Ci b distribution of the flapped portion of the wing. 
The average lift coefficient was estimated from column 
15 of table II. 

The lift-curve slope was calculated from equation 
(3), the value of / being taken from figure 3: 



a=0.999 



0.091 



14 



57.3X0.091 



^0.071 



tX6 



The value of a Q (equation (4)), values of a 0 and a 0f 
having been taken from figures 18 and 19, is 

«o=f|^X0.085+(l-^^X0.099=0.091 
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Figure 19.— Variation of section data with thickness. The N. A. 0. A. 230 series 
airfoils with 0.2c plain flap down 20°; effective Reynolds Number, 8,200,000. 
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Aerodynamic center and pitching moment. — The 

x position was found directly from equation (6): 

^=HA tan A=0.214X6X0.1703=0.219 

The appropriate values of A and H (fig. 5) were used 
and the aerodynamic centers of the wing sections were 
considered to be at the quarter-chord points. 

The pitching-moment coefficient C ms due to the 
moments of the airfoil sections (equation (9)) was 
obtained, as shown in figure 20, from the area above the 
c ma c c 2 curve. Figure 20 illustrates the general method 
that may be applied to any plan form and any distri- 
bution of c ma c across the span. The values of the 
necessary pitching-moment coefficients were taken from 
figures 18 and 19 and are for a value of Ci of zero, in 
which case c m<t c — c m{a c v If c m(a varies appreciably 

with c h it should be taken at the average value of c { over 
the portion of the wing with flaps. In this example, C ms 
coidd also have been calculated from the values of E 
and E' given in figure 7 because the increment, as well 
as the initial value of the pitching-moment coefficient, 
was substantially constant across both the flapped and 
the unflapped parts of the span. 

The pitching-moment coefficient due to the basic lift 
distribution is given by formula (7) 

C mi =GAdA tan A=0.029X0.77X6X0.1703=0.023 

The value of Q was taken from figure 6 and the value of 

Ac i was taken at an intermediate c% (ci=1.0) from 

figure 10 for the N. A. C. A. 23012 airfoil with flap. 

Although Ac i varies with the c t at which it is taken, 

the exact value used does not affect the value of C mi 

h 

appreciably unless the sweepback is large. When the 
quarter-chord points do not lie on a straight line so that 
the angle of sweepback cannot be measured, C mib may 
be computed from formula (8).. The total pitching- 
moment coefficient about the axis through the aero- 
dynamic center is then 

C^ e =C ma +C mi 6 --0.128+0.023 = -0.105 

Drag.— The induced drag was calculated from for- 
mula (13) using values of u, v, and w from figure 8 and 
a value of Acj at an intermediate value of c x (c { =1.0) 
for the N. A. C. A. 23012 airfoil with flap (fig. 10). 
Thus G D t for any value of C L is 



> t Tl> ^gg+(-0.0010)0.77(^-l-0.0100(0.77) 2 

Values of C Di were calculated for a series of values of C L . 

The profile-drag coefficient of the wings was calculated 
by an integration of the section profile drag along the 
semispan as given by 



This integration has been graphically performed as 
shown in figure 21. The value of c do at any point will, 
of course, depend upon the airfoil section, the lift 
coefficient, and the Reynolds Number at that point. 
The calculations are illustrated in table II for a C L of 
0.8 and follow the method of reference 4. 
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Root Distance along semispan, 
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Figure 20— Graphical determination of Cm s for the N. A. C. A. 23012 airfoil with 
0.56 plain flap. 
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Figure 21.— Graphical determination of Cd 0 for the N. A. C. A. 23012 tapered airfoil 
with 0.56 plain flap. 

C D) -|J; C d 0 cd(^)-0.0184 

Values of c dQ c were calculated at intervals along the 
semispan using the known lift distributions. The 
values of (c w ), w , (c* 0|B<n ) m > c h v v and Ac w that 
appear in the various columns of table II were obtained 
from figures 18, 19, and 22. Section data for other 
flap deflections may be found in references 6 and 7. 
The, values of c d0m . n we re extrapolated to the Reynolds 

Number of each point along the semispan by the method 
given in figure 23 of reference 4. Although the formula 
given in reference 4 was derived for sections with 
moderate camber, it should apply approximately to 
airfoils with flaps. 

The Ci distribution for O&=0.8 was obtained from the 
equation 

Ci=C L c lal +c tl) 

where 



C ' al = ~ h L * and C h z 



cb 
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4 5 6 10,000,000 



2 3 4 5 6 10,000,000 6 
Reynolds Number 

Figure 22.— Scale-effect corrections for Ci moi . In order to obtain the section maximum lift coefficient at the desired Reynolds Number, apply to the standard-test value 
the increment indicated by the curve that corresponds to the scale-effect designation (types B, O, D, or E) of the airfoil. (See reference 6, p. 32 and table II.) 

The values of 3* and X 6 at &,/&— 0.5 were obtained from 
reference 2 by cross-plotting against b f /b. The value of 
Ac i at Ci=1.0 for the flapped section was obtained from 
figure 10. 

The profile-drag coefficient is given by the equation 

Cd ° =C *° m u +ACd o 
where Ac dQ depends on the quantity in column 19. of 
table II. Values of Ac ao were obtained from figure 23 
for the sections without flaps (data taken from reference 
6) and from figure 24 for the N. A. C. A. 23012 section 
with flap. 




Figure 23.— Generalized variation of &u a for airfoil sections with flaps neutral. 
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The values of c dQ c given in column 22 of table II 
were plotted in figure 21 and extrapolated to the. flap 
end. ' From the area under the curve, C Dq was found to 
be 0.0184. The process was repeated for other lift 
coefficients and for other wing-and-flap combinations 
and the results are plotted in figures 12 to 17 as effective 
profile-drag coefficients 
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Figure 25.— Calculation of the Cl at which the N. A. C. A. 23012 airfoil with 0.56 
plain flap begins to stall. 



and the value of C Lmax as indicated on figure 25. The 
dashed curve faired through the flap end was used. 
The solid line, which passes through calculated values of 
Ci max —Ci b , would have indicated stalling at the plain 
section just outboard of the flap end at a low C Lmax . 
Observations of the action of tufts, however, indicate 
that stalling does not necessarily begin at this point. 
It appears to be preferable to fair C\ max ~ Ci b through the 
flap end, as shown. The calculated C Lmax value is then 
higher and in better agreement with the test value. 
Tuft observations of the wings with 0.36 and 0.76 
flaps indicated that stalling began at a point other than 
the predicted point so that the method can be expected 
to give only a rough indication of C Lmax . 

COMPARISON OF THE CALCULATED AND THE 
EXPERIMENTAL RESULTS 

The calculated and the experimental results are 
compared in table I. The angles of zero lift and the 
lift-curve slopes are in good agreement. The x positions 
of the aerodynamic center are in fair agreement although 
the experimental aerodynamic-center positions move 
more and more ahead of the calculated position as the 
flap length is increased, probably owing to the forward 
movement of the position of the aerodynamic center 




Flap span /wing span, bf/b 
Figure 26.— Increase in induced drag due to the addition of flaps of various lengths. Cl, 1 .0; Ad, 1 .0. 



Maximum lift coefficient. — The lift coefficients at 
which the wings should begin to stall were estimated by 
the method used in reference 1 except that, instead of 
plotting Ci for a series of values of 0 L , the point at which 
which a section lift coefficient reaches c tmax (c t curve 
becomes tangent to the c lfnax curve) was found more 
conveniently by first deducting c lb from c tmax as in 
figure 25. (See column 23, table II.) The point at 
which the c Zai curve (c ia for (7/ = 1.0) would become 
tangent if expanded to other lift coefficients then deter- 
mines the point where stalling is predicted to start. 

This point is most easily found by calculating 



at several points along the semispan. The minimum 
value gives the location of the predicted stalling point 



of the sections with small flap deflections. The pitching- 
moment coefficients are in good agreement except for 
the N. A. C. A. 5-10-16 wing with 0.56 flap. 

The C De curves given in figures 13 to 17 are in best 
agreement in the region of G Dp . The divergence 

° e min 

for higher and lower lift coefficients is more for these 
wings than for wings without flaps (reference 4). 
Two values of 0 D/> are listed in table I, C Dp and 

C Dg at C L —0.7. It is interesting to note that, for the 
N. A. C. A. 23012 wing, C De increases with flap length 
up to bfjb =0.5 but is then substantially the same at 
b f /b=0.7 as it is at 0.5. The reason for this variation 
is that the increase in profile drag with flap length is 
compensated by the reduction in induced drag beyond 
b f fb=0.5. If plain flaps at a moderate angle are used 
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Wing lift coefficient, C L 
Figure 27.— Percentage change in induced drag relative to that for flaps neutral. 0.56 plain flap; A, G. 



Percentage change given by 



CL 2 /rAu 



X100. 



for take-off, as long a flap span as possible is therefore 
indicated to obtain the lowest drag. The decrease in 
induced drag with flap length beyond b f /b equal to 
about 0.36 is illustrated in figure 26. The curves are 
given for various taper ratios, two values of A, and for 
C L and Act both equal to 1.0. A value of Ac { =1.0 
corresponds to an airfoil of moderate thickness with a 
0.2c flap at an angle a little larger than 20°. The 
change of induced drag with lift coefficient is illustrated 
in figure 27 for a wing with half-span flap of aspect 
ratio 6 and for two values of Acj. 

The calculated and the experimental C Lopt values 
given in table I are in as good agreement as can be 
expected in view of the difficulty of determining 
c lgpt . The agreement of the C Lmas values is fair. 



CONCLUDING REMARKS 

Although test results and the comparison with 
calculated results have been given only for the case of a 
plain flap deflected 20°, other test results were available 
for the N. A. C. A. 23012 wing with a split flap 45° 
down and with the same flap spans as those used herein. 
Comparison of the test and the calculated results for 
the split flap showed agreement similar to that obtained 
for the plain flap deflected 20°. It thus appears that 
fair estimates of the characteristics of tapered wings 
with partial-span flaps deflected various amounts can 
be obtained from the factors and the method given. 



Langley Memorial Aeronautical Lap oratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., January 23, 1939. 



APPENDIX A 



SYMBOLS USED IN TEXT 

Cj, section lift coefficient. 
Ac h increment of section lift coefficient due to 

flap deflection. 
Ci b , section basic lift coefficient (C L =0). 
Ci a , section additional lift coefficient. 
€i al , section additional lift coefficient for C L = 1.0. 
Ca 0 , section profile-drag coefficient. 
(i=o)' wm S angle of attack for zero lift, measured 

from the root chord. 
ai Q , angle of zero lift of root section. 

b, wing span. 

b f} total flap span. 
S, area of wing. 

S f , area of part of wing equipped with flaps. 
A, wing aspect ratio, b 2 /S. 
X, taper ratio, c t /c s . 

q, dynamic pressure. 

c, chord at any section along the span. 

c t , tip chord (for rounded tips, c t is the fictitious 
chord obtained by extending the leading 
and the trailing edges to the extreme 
tip). 

c s , chord at root of wing or plane of symmetry. 
A, angle of sweepback measured between the 
lateral axis and a line through the aero- 
dynamic centers (approximately the 
quarter-chord points) of the wing sec- 
tions. 
S T , flap angle. 
C L , wing lift coefficient. 
C D , wing drag coefficient. 
Cz) 0 , wing profile-drag coefficient. 
C Df , effective wing profile-drag coefficient. 
C D{ , wing induced-drag coefficient. 

16 



c mae , section pitching-moment coefficient about 
section aerodynamic center. 
c m e )0 , section pitching-moment coefficient about 
aerodynamic-center position with flap 
, neutral. 

c m „, section pitching-moment coefficient with 

flaps neutral. 
Ac m , increase in section pitching-moment coeffi- 
cient above c mQ due to flap deflection. 
M, total wing pitching moment. 
Mi b , wing pitching moment due to basic-lift 
forces. 

C mi , wing pitching-moment coefficient due to 
b 

basic-lift forces. 
C ms , wing pitching-moment coefficient due to the 
pitching moments of the wing sections. 
@m a c , total wing pitching-moment coefficient about 
aerodynamic center. 
a, wing lift-curve slope. 
a 0 , lift-curve slope of section without flap. 
a 0p lift-curve slope of section with flap. 
x, moment arm measured from the quarter- 
chord point of the root chord and parallel 
to it (positive rearward). 
y, lateral distance. 

y f , lateral distance to inboard end of flap. 
Xa.c, coordinate of wing aerodynamic center. 
R e , effective Reynolds Number. 
L a , additional load parameter. 
L b , basic load parameter. 
J, factor of angle of zero lift. 
H, factor of wing aerodynamic center. 
G, factor of basic-lift pitching moment, 
jf, factor of wing lift-curve slope. 
E and E' f factors of section pitching moment. 
u, v, w, factors of induced drag. 



APPENDIX B 



AERODYNAMIC FACTORS IN TERMS OF THE FOURIER COEFFICIENTS 

The various aerodynamic factors were obtained from 
a Fourier analysis in which the circulation r was 
expressed (see reference 2) by 

z » 

where 

c s is the chord at plane of symmetry. 
m„ slope of the section lift curve at the plane of 
symmetry, per radian. 
V, wind velocity. 

cos «==f- 

If the Fourier coefficients of the plain wing at an 
angle of attack of one radian are denoted by A n and if 
the Fourier coefficients for the same wings with a 
constant angle of attack extending only over the 
center of the span are denoted by a n , the various 
aerodynamic factors (in terms of the Fourier coefficients) 
can be found from the following equations: 

Aim Q 

in which m 0 is the slope of the lift curve at any section, 
per radian. 

rr__ 2 / Ai A3 fls, Ar Ag \ 

5 21 i "45 77 ' ' V 

r _c*AT(<h <h , <h 1 \ <h_/A ? _As,A ! . \~| 
26L\5 2l" t "45 i " ' V A l \ 5 21 i ~45 i ~* " /J 
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TABLE I 

COMPARISON OF CALCULATED AND EXPERIMENTAL CHARACTERISTICS 
[Wings with 0.2c plain flaps 20° down) 



Wing 


Flap 
length 
(frac- 
tion 
span) 


X= 

Cl/C, 


Aspect 
ratio 


Boot sec- 
tion, 
N. A. C. A. 


Construc- 
tion tip 
section 
(extreme 
tip), 
N\ A. C. A. 


Sweep- 
hack 
(deg.) 


a, (L-0> 

(deg.) 


a 


x a . t . 
S/b 

(0 


Experi- 
mental 


Calcu- 
lated 


Experi- 
mental 


Calcu- 
lated 


Experi- 
mental 


Calcu- 
lated 


Tapered N. A. C. A. 23012.. 
Tapered N. A. C. A. 5-10-16. 


I'i 

{'.. 


J 0.5 

I •» 


6 

10 


23012 
23016 


23012 
23009 


9.67 
0 


( -1.3 
-4.8 
1 -6.9 
I -8.7 
f -1.2 
\ -6.2 


-1.2 
-4.6 
-6.7 
-8.6 
-1.1 
-6.4 


0.075 
.072 
.070 
.067 
.083 
.079 


0.076 
.073 
.071 
.069 
.083 
.078 


0.210 
.209 
.201 
.193 
-.011 
-.021 


0.219 
.219 
.219 
.219 

0 
0 


Wing 




Co. 




CL mat 


Cd, at 


Cl=0.7 


Experi- 
mental 


Calculated 


Experi- 
mental 


Calculated 


Experi- 
mental 


Calculated 


Experi- 
mental 


Calculated 


Experi- 
mental 


Calculated 


Tapered N. A. C. A. 23012— 
Tapered ST. A. C. A. 5-10-16. 


f -0.014 
-.068 
1 -. 106 
I -. 143 
f -.009 
\ -.136 


-0.008 
-.068 
-.105 
-.142 
-.007 
-.157 


0. 0076 
.0165 
.0175 
.0156 
.0080 
.0190 


0.0071 
.0155 
.0170 
.0158 
.0081 
.0178 


0.13 
.20 
.20 
.16 
.03 

-.12 


0.09 
-.02 
.03 
.05 
.10 
-.23 


1.71 
1.83 
1.91 
1.98 
1.50 
1.77 


1.67 
1.77 
1.81 
1.90 
1.49 
1.74 


0. 0106 
.0182 
.0205 
.0203 
.0123 
.0262 


0. 0100 
.0209 
.0226. 
.0222 
.0113 
.0274 



i Measured from the quarter-chord point of the root chord, positive toward the trailing edge. 

Table II 

CALCULATION OF C Do FOR C'z,=0.8 
[N. A. C. A. 23012 tapered wing with 0.56 flap deflected downward 20°; R e (based on S/b) =8,200,000] 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


Distance 
from 
center, 
fraction 
semispan, 
V 

6/2 


thickness 
chord 

t/c 


Chord 
c 

(in.) 


Effective 
Reynolds 
Number 

R. 
(millions) 


at 

i2, =8,200,000 


Ac, »<ij 

(fig. 22) 


(5) +(6) 


("0 M ,n),« d 

at 

R.= 8,200,000 




C 'opl 




Load 
parameter 

La 

(reference 2) 


0.00 
.15 
.30 
.45 
.60 
.70 
.80 
.90 
.95 


0.12 
.12 
. 12 
.12 
.12 
.12 
.12 
.12 
.12 


6.82 
6.31 
5.80 
5.28 
4.77 
4. 43 
4. 09 
3.50 
2.66 


11.18 
10.35 
9.51 
8.60 
7.82 
7.26 
6.71 
5.74 
4.36 


2. 06 
2. 06 
2.06 
2.06 
1.75 
1.75 
1.75 
1.75 
1.75 


0.02 
.02 
.01 

0 

-.01 
-.02 
-.03 
-.06 
-.10 


2.08 
2.08 
2.07 
2.06 
1.74 
1. 73 
1.72 
1.69 
1.65 


0. 0128 
.0128 
.0128 
.0128 
.0071 
.0071 
.0071 
.0071 
.0071 


0.0125 
.0126 
.0127 
.0128 
.0072 
.0072 
.0072 
.0074 
.0076 


0.10 
.10 
.10 
.10 
.10 
.10 
.10 
.10 
.10 


1. 98 
1.98 
1.97 
1.96 
1.64 
1.63 
1.62 
1.59 
1.55 


1.291 
1.263 
1. 191 
1.107 
.995 
.908 
.789 
.607 
.447 



1 


13 


14 


15 


16 


17 


18 


19 


20 


21 


22 


23 


y 

6/2 




Load 
parameter 
Lt 

(reference 2) 


ACiS 


CiX.ci ai = ci a 


CI 


C '~ Cl opt 






Cd 

(9) +(20) 


Cd 0 C 








0.00 
.15 
.30 
.45 
.60 
.70 
.80 
.90 
.95 


0.946 
1.001 
1.027 
1.048 
1.043 
1.025 
.964 
.867 
.840 


0.289 
.276 
.235 
.030 
-.190 
-.252 
-.266 
-.224 
-.168 


0. 163 
.168 
.156 
.022 
-. 153- 
-.219 
-.250 
-.246 
-.243 


0. 757 
.801 
.822 
.838 
.834 
.820 
.771 
.694 
.672 


0. 920 
.969 
.978 
.860 
.681 
.601 
.521 
.448 
.429 


0. 820 
.869 
.878 
.760 
.581 
.501 
.421 
.348 
.329 


0.414 
.439 
.446 
.388 
.354 
.307 
.260 
.219 
.212 


0.0U7 
.0127 
.0129 
.0106 
.0023 
.0018 
.0013 
.0010 
.0009 


0.0242 
.0253 
. 0256 
.0234 
.0095 
.0090 
.0085 
.0084 
.0085 


0. 1650 
. 1596 
.1485 
.1236 
.0453 
.0399 
.0348 
.0294 
.0226 


1.92 
1.91 
1.91 
2. 04 
1.89 
1.95 
1.97 
1.94 
1.89 
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Positive directions of axes and angles (forces and moments) are shown by arrows 



. : — "T* — '— 

Axis 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designation 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal 


X 
Y 
Z 


X 
Y 
Z 


Rolling. 
Pitching.... 
Yawing- _ 


L 

M 

N 


Y +Z 

Z >X 

x — >r 


Roll 


4> 
6 

4> 


u 

V . 
w 


V 
<l 
r 


Lateral .__ 

Normal- _ 


Pitch..,_ 
Yaw. __ 







Absolute coefficients of moment 

n-JL r -Ml 

° l ~qbS W gcS 

(rolling) (pitching) 



D y Diameter 
P, Geometric pitch 
p/D, Pitch ratio 
V, Inflow velocity 
V s , Slipstream velocity- 
Thrust, absolute coefficient C T 



N_ 
gbS 
(yawing) 



Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 



4. PROPELLER SYMBOLS 



T, 



T 

Q, Torque, absolute coefficient C Q — 9ps 



P, 

v, 
n, 



j 

Power, absolute coefficient Cj>= 3n § 

Speed-power coefficient — -W p-§ 
Efficiency 

Revolutions per second, r.p.s. 
Effective helix angle— tan~ J jfs-^— ^ 



1 hp.=76.04 kg-m/s=550 ft-lb./sec. 
1 metric horsepower =1.0132 hp. 
1 m.p.h. =0.4470 m.p.s. 
1 m.p.s. =2.2369 m.p.h. 



5. NUMERICAL RELATIONS 
1 lb 



0.4536 kg. 
1 kg=2.2046 lb. 
1 mi. = 1,609.35 m= 5,280 ft. 
1 m=3.2808 ft. 



